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FOREWORD 


This document is based on a thesis submitted by the author to the faculty of 
the George Washington University School of Engineering and Applied Science. 
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CLOSE- LOOP DOPPLER CORRECTION SCHEME 


Peter E. O'Neill, Jr. 
Communications Research Branch 


ABSTRACT 

Doppler frequency shift occurs when there is relative velocity be- 
tween a satellite and ground station. The normalized frequency error 
is v /c , where v is the relative velocity between the satellite and the 
ground station and c is the speed of light. 

This document describes a technique for reducing the effective 
Doppler frequency shift. This technique can reduce the Doppler effect 
to approximately second order, (v/c ) 2 . Experimental data have verified 
the technique, differing from theoretical predictions by 3 percent. A 
ground-controlled satellite clock with a stability of 5 parts in 10 10 was 
established in Relay n. Without the Doppler correction scheme, the 
normalized frequency error can exceed 1 part in 10 5 . 

This document presents what is believed to be a unique analysis of 
the Doppler correction scheme in the variable velocity case with as- 
sociated experimental verification. 
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CLOSED- LOOP DOPPLER CORRECTION SCHEME 


INTRODUCTION 

When there is relative motion between a transmitter (T) and a receiver (R), 
the frequency of the signal detected at the receiver will not be the same' as that 
sent from the transmitter. This change in frequency is called the Doppler fre- 
quency shift. 

By using closed-loop techniques, it is possible to correct for Doppler fre- 
quency shift. The signal at R is retransmitted back to T . The frequency dif- 
ference between this signal and the originally transmitted signal is a measure of 
the two-way Doppler frequency shift. Using the standard technique , 1 it can be 
shown that: 


f 


R 





c +v 


l R’ 


c 

c +v 



C - V 

c +v 


f 


T * 


where f is the frequency received at R, 

f T is the frequency transmitted at T , 
f ' is the frequency received at T. 
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Brown, R. C.: A Textbook of Physics. Longmans, London, 1961, p. 1013 
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Because v « c , 


fl ^ c -2v 

*R ~ 

K C 


(4) 


For small values of v /c , the frequency shift of the signal received at R will be 
approximately one-half the frequency shift of the signal received at T . 

A servomechanism was devised which forced 

f R + f T = 2f REF , (5) 

where f REp is the desired signal frequency in the absence of Doppler. Substituting 
equation (4) into equation (5) and solving for f T yields: 


c - 2v 


f T + f T 2 f REF , 


f = f 

T r REF 


( 6 ) 


By using equation (1), it is now possible to determine the frequency which will 
be received at R. 


f = c _v f 

K “ ” r T ’ 


C -V C 


C C - V 


f 

REF ’ 


f = f 

R REF 


(7) 
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Appendix A contains a description and block diagram of the servomechanism. 

It has been shown that the Doppler frequency shift can be eliminated from a 
signal under the following special circumstances: 


1. v « c 

2. v t v (t) 

If v = v(t), then the two v’s in equation (3) are not equal, because they are evaluated 
at different times. 

The preceding development is common in the literature on this subject. The 
case where v = v(t) has not been thoroughly analyzed. The case where v = v(t) 
is treated both theoretically and experimentally in this document. 


ANALYSIS OF THE CONSTANT VELOCITY CASE 

In order to keep the analysis as general as possible, the constant velocity 
case will be rederived without the constraint that v « c . Repeating equation (1) 
and equation (3), 


( 1 ) 
( 3 ) 

where f R is the frequency received at the satellite, 

• • 

f T is the frequency transmitted from the ground, 
fg is the frequency received on the ground. 

The servomechanism forced 



c . t - if 
*R T *T “ ‘ REF ' 


( 8 ) 
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where f REF is the desired signal frequency in the absence of Doppler. Substituting 
equation (3) into equation (8) yields: 


c - v 
c +v 


f- = 2f 


REF 1 


( 9 ) 


f 


T 


C +V 
C 


REF' 


( 10 ) 


By substituting this into equation (1), it is possible to find the signal which will be 
received at the satellite. 



( 1 ) 


C - V C + V 


1 REF ’ 



( 11 ) 


Note that the Doppler is not reduced to zero as before, but is reduced to 
second order in v /c . 

A normalized frequency error is defined for use in subsequent sections. 


En = 


_ f 

REF 

f 

REF 


( 12 ) 


= - (v/c) 2 . 


(13) 
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ANALYSIS OF THE VARIABLE VELOCITY CASE 
Classical Approach 

This section describes a still more general analysis of the closed-loop Doppler 
correction scheme in which v/c is significant and velocity is a function of 
time. Let 


v = v(t) , 


= V + at , 


(14) 


where v is the velocity at time t = 0, 

a is the acceleration at time t = 0. 

For short periods of time (t < 100 sec), the relative velocity between a ground 

station and satellite can generally be expressed as the summation of an initial 
velocity and a constant acceleration. 

The frequency received at the satellite at time t , f R (t) , is a function of the 
relative velocity at time t , and the frequency transmitted at time t-T, f T (t -T). 
T is the finite propagation time of the signal from the ground to the satellite, and 
is also a function of time. Substituting equation (14) into equation (1) yields: 


*«<«> 


t- v - « f T (t--r). 

c 1 


(15) 


The frequency received at the ground at time t , f R (t) , is a function of the 
relative velocity at time t-T and the frequency transmitted from the satellite 
at time t - T , f R (t - T). T is the finite propagation delay and is equal to T in 
equation (15) because the turnaround time of the signal in the satellite is negligible 
relative to the propagation time delay. Substituting equation (14) into equation (2) 
yields: 


f R (*> = 


° 7 T “ t n M 1 - 1 )* 

c +v + a (t -T) 


/-* \ 
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Replacing t with t - T in equation (15) and substituting into equation (16) 
yields: 


f R(*> 


c - v - a (t -T) 
c + v + a (t -T) 


f T (t-2T). 


( 17 ) 


Equation (17) provides a relationship between the frequency received on the ground 
at time t , (t) , and the frequency transmitted from the ground at time t - 2T, 

f T (t - 2T) . To describe the servomechanism, however, it is necessary to have 
a relationship between f T (t) and f R (t) . Because there was no obvious way to 
establish this relationship, this approach was abandoned. 


Series Approach 


The equations of motion of most macroscopic objects (e.g., a satellite) gen- 
erate relatively smooth, continuous curves. Short segments of smooth continuous 
curves can be accurately represented by a truncated power series. The range 
to a satellite, for example, can be defined as: 


R (t) = r + vt + I-at 2 . 


(18) 


where 


R (t) is the range at any time, 
r is the range at time t = 0, 
v is the velocity at time t = 0, 
a is the acceleration at time t = 0. 


This equation is valid only over a finite time interval because the satellite is 
still at a finite range after infinite time. Over a proper time interval, however, 
the errors introduced by the use of equation (18) are insignificant. 

Because the Doppler effect is linearly related to the rate of change of range, 
the Doppler frequency perturbation can also be adequately represented by a 
truncated power series. The phase at the satellite, at any time t , can be ex- 
pressed as: 


<£(t) = b Q + (bi + w) t + b 2 t 2 , 


(19) 
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where 


4 > (t) is the phase at any time t , 

b Q is the initial phase at time t = 0, 

(b x + w) is the initial frequency at time t = 0, 

b 2 is the initial frequency rate at time t = 0. 


In equation (19), b x is the frequency error in radians per second caused by 
Doppler, and w is the frequency in radians per second in the absence of Doppler 
(2tt f REF ). For simplicity, assume that b 2 is small relative to b x + w (this 
assumption is justified by the experimental data) and therefore has negligible 
effect on0(t) for small values of t. Hence, 


<£(t) = b o + ( b x + w) t , (20) 

where b x is a function of both the initial velocity and acceleration. 

The phase of a signal being transmitted from the ground at time t is the 
identical phase which will be received at the satellite at time t + T, if T is the 
propagation time from the ground to the satellite of a signal arriving at the 
satellite at time t + T. Therefore, 


^(t) = b 0 + ( b x + w) (t + T), 


( 21 ) 


where 0 T (t) is the phase of the signal transmitted from the ground at any time 
t in terms of the phase at the satellite. 

The phase of a signal being received on the ground at time t is the identical 
phase transmitted from the satellite at time t - r , if r is the propagation time 
from the satellite to the ground of a signal transmitted from the satellite at 
time t - r . Therefore, 


<£r (t) = b 0 + ( b i + w) (t - r), 


( 22 ) 


where 0 R (t) is the phase of the signal received on the ground at time t in terms 
of the phase at the satellite. 
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The servomechanism causes the sum of 4> T (t) and ^ (t) to equal a reference, 
or 

(t) + 0R (t) = 2wt + (23) 

where <// is an arbitrary initial phase angle. Substituting equation (21) and equation 
(22) into equation (23) yields: 

b o + ( b i +w ) (* +T ) + b 0 + ( b i +w ) (t -t) = 2wt +\p, (24) 

2b Q + 2 (bj +w)t + (b t + w) (T- t) = 2wt + \p. (25) 

Both T and r can be expanded into a Maclaurin series so that 

T(t) = T(0) + T(0)t +T(0) |i, (26) 

r(t) = r(0) +f(0)t +r (0)l!. (27) 

Substituting this into equation (25) yields: 

2b 0 + 2 (b x +w)t + (b 1 +w) T (0) + T (0) t 

+ T (O)li-r(O) -f(0)t - V (0)|i 

= 2wt + (28) 
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Equating powers of t yields: 


2b 0 + (b, +w) [T (0) — t (0) ] = 0 , 

(29) 

(bj +w) [2 + T(0) - r (0)] = 2w, 

(30) 

(bj +w) [T (0) - r(0)] = 0. 

(31) 


Solving equation (30) for yields 


. -w [T CO) - r (0)] 
2 + T (0) - f (0) 


(32) 


where is the frequency error in radians per second due to Doppler in the 
variable velocity case. The general expressions for T (t) and r (t) , as well 
as their series expansions, are developed in Appendix B. Substituting the ex- 
pressions for T (t) and f (t) into equation (32) yields: 


bj = -w 


(c-v)v + ra (c+v)v-ra 

(c - v ) 2 (c +v ) 2 

2 + (c-v)v + ra (c +v)v - ra 

(c - v ) 2 (c + v ) 2 _ 


(33) 


In Appendix C, a first-order approximation is made for equation (33): 


b x ~ W (v 2 + ra). 
c 2 


(34) 
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The normalized theoretical error, E NT , is obtained by dividing equation (34) by w. 


_ -(v 2 + ra) 


"'NT 


(35) 


This is in complete agreement with the constant velocity case. If the acceleration, 
a , in equation (35) is set equal to zero, then the normalized error calculated by 
using equation (35) will equal the normalized error found by using equation (13). 


EXPERIMENTAL TEST SYSTEM 

General Description 

A comparison of the corrected frequency at the satellite with the reference 
frequency was required to experimentally determine the quality of the Doppler 
correction scheme. However, the satellite was only a simple repeater and con- 
tained no reference clocks. It was therefore impossible to make the required 
measurements in the satellite. 

To overcome this difficulty, the system was modified so that the signal 
looped through the satellite twice. This placed the center of the system, a 
virtual satellite, on the ground. A description and block diagram of the servo- 
mechanism of the test system is contained in Appendix D. 

Because both the corrected signal at the virtual satellite and the reference 
frequency were available on the ground, it was possible to measure the quality 
of the Doppler correction scheme. 

Appendix D also contains a description and functional block diagram of the 
method used to compare the frequency at the virtual satellite, f vs , with the 
reference frequency, f pgTi . . 


Test System Analysis 


The analysis of the test system is similar to that of the operational system. 
Assume that, as in equation (20), the phase at the satellite, 0 (t) , can be expressed 
as 


= b 0 + ( b i + w ) t > 


(36) 
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where 


4>(t) is the phase at any time t , 

b 0 is the initial phase at time t = 0, 

(b x +w) is the initial frequency at time t = 0. 

As defined in the section on the analysis of the variable velocity case, b t is the 
frequency error in radians per second caused by Doppler and w is the frequency 
in radians per second in the absence of Doppler (2tt f REF ). 

The transit time of a signal from the ground to the satellite is exactly the 
transit time of the same signal from the satellite to the ground if it is retrans- 
mitted immediately. This is intuitively justified by noting that the range up is 
equal to the range down if the signal is retransmitted immediately. The phase 
of a signal being transmitted from the ground at time t is the identical phase 
which will be received at the virtual satellite at time t + 2T , if T is the propaga- 
tion time from the ground to the actual satellite of a signal arriving at the actual 
satellite at time t + T. Therefore, 


<Mt) =b 0 + ( b i+ w ) (t+2T), 


(37) 


where 4> y (t) is the phase of the signal transmitted from the ground at any time 
t in terms of the phase at the virtual satellite. 

The phase of a signal being received on the ground at time t is the identical 
phase which was transmitted from the virtual satellite at time t - 2 t, if r is the 
propagation time from the actual satellite to the ground of a signal transmitted 
from the actual satellite at time t - r . Therefore, 


=b 0 + ( b l +W > (* - 2T )> 


(38) 


where 0 R (t) is the phase of the signal received on the ground at time t in terms 
of the phase at the virtual satellite. 

The servomechanism causes the sum of 4> y (t) and <££ (t) to equal a reference, 
or 


<t> T (t) + (t) = 2wt +0, (3y) 

where is the arbitrary initial phase angle. Substituting equation (37) and 
equation (38) into equation (39) yields: 
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t»o + (bj +w) (t + 2T) + b 0 + (b t +W) (t - 2 r) = 2wt + \p , 


(40) 


2b Q + 2 (bj +w) t +2 (bj +w) (T - r) = 2wt + <//. 

(41) 

As in the preceding section, T and r can be expanded into a Maclaurin series so 
that 

T (t) = T (0) + T (0) t + T (0) t l, 

(42) 

T(t) =T( 0) + t (0)t + V(0) ii. 

(43) 

Substituting this into equation (41), 


2b 0 + 2 (bj +w) t + 2 (bj +w) j^T (0) + T (0) t 


+ T (O)li-r(O) -r(0)t - r (0) ^.] 


= 2wt + i/». 

(44) 

Equating powers of t yields: 


2b 0 + 2 (b t +w) [T (0) - r (0)] = if* , 

(45) 

(bj +w) [1 + T (0) - r (0)] = w. 

(46) 

(bj + w) [T (0) - r(0>] = 0. 

(47) 
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Solving equation (46) for b t yields : 


b 


l 


-w[T (0)-f(0)] 

. .. n . 9 

1 + T (0) - f (0) 


(48) 


where b x is the frequency error in radians per second due to Doppler. Appendix 
B contains the general expressions for T (0) and r (0) . Substitution yields: 


bj = -w 


(c -v)v + ra (c + v) v - r a 

(c - v) 1 2 3 (c + v) 2 

1 + (c - v) v + ra (c + v) v - r a 
(c - v) 2 (c +v) 2 


(49) 


In Appendix E, a first-order approximation is made for equation (49) •' 


-2w (v 2 + r a) 


(50) 


The normalized theoretical error of the test system, E' NT , is obtained by dividing 
equation (50) by w . 


e: 


NT 


-2 (v 2 + ra) 
^2 


(51) 


A comparison of equation (51) and equation (35) indicates that the error at the 
virtual satellite in the test system will be exactly twice the error at the actual 
satellite in the operational system. 


1 till X 

The Relay satellite was used to obtain experimental data. Figures 1, 2, and 

3 show the range, velocity, and acceleration of Relay with respect to the Mojave 
ground station as a function of time after acquisition for orbit 4472, September 17, 
1965. 
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The first experimental run was initiated at time t = 2.75 minutes with the 
following initial conditions : 

r = 6.8 x 10 6 meters, 
v = 4.2 x 10 3 meters per second, 
a = - .25 meters per second squared. 

The second experiment was started at time t = 25.25 minutes with the 
following initial conditions : 

r = 3.3 x 10 6 meters, 
v = 2.0 x 10 3 meters per second, 
a = 4.3 meters per second squared. 

Note that v 2 » ra for the first part of the experiment and ra > v 2 for the 
second part, thus providing a reliable test of the theoretical prediction. 

The recorded data are in the general form of a scatter diagram as shown in 
Figure 4. A least-squares polynomial approximation is made for the data as 
discussed in Appendix F. 2 

Table 1 lists the normalized measured errors, E NM , as determined by the 
computer and the normalized theoretical errors for the test system, E' NT , as 
calculated in Appendix G. Because of the close agreement between the measured 
and predicted results, it is assumed that the theoretical analysis is valid. 


Table 1 

Comparison of E NM * and E' NT ,„ 





Difference 

Test No. 

®NM 

E'nt 

Between E NM and E' NT 
<%> 

1 

-3.59 parts in 10 10 

-3.54 parts in 10 10 

1.4 

2 

-7.40 parts in 10 10 

-7.16 parts in 10 10 

3.0 


*Normalized measured error of the test system. 
**Normalized theoretical error of the test system. 


2 

Hilderbrand, F. B.: Introduction to Numerical Analysis, ch. vii, McGraw-Hill Book Co., Inc., 
New York, 1965 
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Figure 1 -Satellite Range vs Time 
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Figure 2-Satellite Velocity vs Time 
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ACCELERATION (METERS PER SECOND SQUARED) 



TIME AFTER SATELLITE ACQUISITION (MINUTES) 
Figure 3-Satellite Acceleration vs Time 
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1 





TIME (t ) 

Figure 4— Scatter Diagram of Recorded Data 


WORST-CASE ANALYSIS FOR RELAY II 

Relay II, launched in 1964, is in an elliptical orbit with an apogee of 7500 km 
and a perigee of 2100 km. The worst-case relative velocity and accelera- 
tion occurs when perigee is approximately overhead. For a simplified worst- 
case analysis, assume that the orbit is circular at the altitude of perigee, 2100 km, 
and that the orbital angular velocity of the satellite is constant and equal to the 
angular velocity at perigee, 7 x 10 ~ 4 rad/sec. 

Let r e (Figure 5) be the radius of the earth, 6371 km. Let d represent 
the distance from the center of the earth to the satellite. If a ground station 
is located at g on the surface of the earth, then the range to the satellite is: 3 


3 

Handbook of Mathematical Tables. 2d ed., Chemical Rubber Publishing Co., Cleveland, 1964, 
p. 570 
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r = (d 2 + r e 2 - 2<lr e cos 6 ) 1/2 . 


( 52 ) 


But, 


0 = *( t), 


( 53 ) 


□ at, 


where a = angular velocity of satellite referenced to the center of the earth at 
perigee. 


Therefore, 


r = (d 2 + r e 2 - 2dr e cos a t) 1/2 , 


(54) 


v = = dr a sin a t (d 2 + r 2 - 2dr cos a t) _1/2 , 

dt e v e e ' 


(55) 


a = -„1 = dr e a 2 cos a t (d 2 + r 2 - 2dr cos a t)* 1 ^ 2 
dt 2 

-d 2 r 2 a 2 sin 2 a t (d 2 + r 2 - 2dr cos a t) _3/2 . (®6) 

e v e e 7 


In the operational system the normalized error was, from equation (35): 

_ -(v 2 + ra) 

^NT 2 

C 

Substituting, 


^NT 


d 2 r e 2 a 2 sin 2 at 
d 2 + rj 2 - 2dr g cos at 


dr a 2 cos a t 


-d 2 r 2 a 2 sin 2 at 

e 

d 2 + r e 2 - 2dr e cos at 


= - (dr a 2 cos a t)/c 2 . 


(57) 
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When 6 = a t =0, E NT = E NT (max), 


E NT (max) = -dr e a 2 /c 2 , (58) 

where d = 7500 km + 6371 km, 

= 13,871 km, 
r g = 6371 km, 
a = 7 x 10~ 4 rad/sec, 
c = 3 x 10 s km/sec, 


E nt (max) = 4.7 x 10' 10 . ( 59 ) 

In the operational system the worst-case frequency error due to Doppler for 
Relay II is 4.7 parts in 10 10 . Without the Doppler correction scheme, the maxi- 
mum normalized Doppler error would be approximately 1.25 parts in 10 s . With 
the Doppler correction scheme, the maximum error will occur directly overhead 
(6 = 0); without it, the maximum Doppler occurs at the maximum range. 


CONCLUSIONS 

It has been demonstrated that the stability of ground-controlled satellite 
clocks can be greatly increased by the use of the Doppler correction scheme. 

It has been shown both theoretically and experimentally that the correction 
scheme is useful in the variable velocity case. The close agreement between 
the theoretical and measured results (Table 1) indicates that the theoretical 
analysis is valid. 

Maximum error with the Doppler correction scheme for Relay II occurs 
when the satellite is directly overhead. Because the relative velocity of the 
satellite becomes zero when the satellite is directly overhead, the maximum 
error is a function of the product of range and acceleration, equation (35). The 
true significance of this variable velocity analysis is now obvious because the 
constant velocity analysis predicts zero error when the satellite is directly 
overhead, at which time the relative velocity is zero, equation (13). 


The worst-case normalized error for Reiay II was 4.7 parts in 10 10 . With- 
out the Doppler correction scheme, this normalized error would be greater than 
1 part in 10 s . 
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APPENDIX A 

SERVOMECHANISM OF THE OPERATIONAL SYSTEM 

A comparison of equation (1) and equation (4) indicates that, for small values 
of v/c , the Doppler frequency shift of the signal received at the satellite is 
approximately one-half the Doppler frequency shift of the signal received on the 
ground. 


f 


R 



( 1 ) 


'c -2v 


) 


( 4 ) 


where f T is the frequency transmitted from the ground, 

f R is the frequency received at the satellite, 
f R is the frequency received on the ground. 

Therefore, the frequency of the signal at the satellite is the average of the fre- 
quencies of the signals transmitted and received on the ground. If the average of 
the frequencies of the signals transmitted and received on the ground, f T and f R , 
respectively, is maintained at f REF , then f REF will always be present at the 
satellite. 

The servomechanism must therefore perform the following function: 


f + f' 

T + R 


= f 


REF 


(Al) 


or 


*T + *R 2 *REF ’ 


(A2) 


£ £ * r\ £ _ n 

1 T r 1 R — **REF ~ u - 


/d o\ 
v* *•*'/ 
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Figure A-l is a simplified block diagram of the servomechanism. The sum, 
f T and f R , is formed in mixer 1 and is passed through a selective filter to 
mixer 2. 2 f REF is also an imput to mixer 2. The difference output of mixer 2 is 

f R + f T " 2f REF* 



Figure A-l —Servomechanism of Operational System 
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This output is passed through the low-pass loop filter, kf (s) , and is applied as 
a correction voltage to the voltage-controlled crystal oscillator (VCXO). The 
output of the VCXO is directly f T . The correction voltage will be zero only when 

+ *T 2f RE F = 0- 

When f R (t) is instantaneously perturbed by Doppler, 

f'(t) + f T (t)^2f REr . (A4) 


A correction voltage is instantaneously developed at the output of the low-pass 
filter, kf (s) , which adjusts the VCXO so that the sum of f R (t) and f T (t) is 
again 2 f R g F • 

Because the average of f R and f T is maintained at f REF , f REF will always 
be present at the satellite. 
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APPENDIX B 

CALCULATION OF TIME DELAYS 


As discussed in the section on the variable velocity case, the range can be 
expressed as: 


R (t) = r + vt + —at 2 , (Bl) 

2 

where R (t) is the range at any time t , 
r is the range at time t = 0, 
v is the velocity at time t = 0, 
a is the acceleration at time t = 0. 

The range at some other time t + T can be expressed as : 

R (t +T) =R(t) T, (B2) 


where 


dR 

— = v + at. 
dt 


(B3) 


Substituting equation (B3) and (Bl) into (B2), 

R (t + T) = r + vt + i. at 2 + vT + at T. (B4) 


The time, T , required for a signal to reach the satellite is equal to the range 
when the signal reaches the satellite, R (t + T), divided by the speed of light, c . 
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Therefore, 


T _ R (t + T) 
c 


R (t + T) = cT. 


Substitution of equation (B6) into equation (B4) yields: 


cT = r + vt + _at 2 + vT + at T. 
2 


Gathering terms in T , 


(c-v-at)T=r +vt + _ at 2 , 


or 


T = 


r + vt + _ at 2 
2 

c - v - at 


It is useful to expand this into a Maclaurin series of the form 

f (t) = f (0) + f (0) t. 


Therefore, 


T (t) = T(0) + T (0) t, 


where 



(B13) 


T (0) = ( c ~ v ) v +ar 
(c -v) 2 


The range at some other time t - r can be expressed as: 


R(t-T)=R(t)-^-T, 


where 


• = v + at . 


Substituting equation (Bl) and (B15) into (B14), 


R(t-r) = r + vt + A-at 2 - vr-at r. 


The time, r , required for a signal to arrive from the satellite is equal to the 
range when the signal was transmitted from the satellite, R (t - r> divided 
by the speed of light, c . 


T _ R (t - t) 


Therefore, 


R (t -T) = c. 


Substitution of equation (B18) into equation (B16), 


c r = r + vt + i-at 2 - vr-atr. 
2 
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Gathering terms in r , 


(c+v+at)r = r+ vt 



(B20) 


or 


r 


r + vt + _at 2 
2 

c + v + at 


Expansion into a Maclaurin series yields: 

T (t ) = T(0) + f (0) t, 

where 


T (0) = -I— , 

c + v 


(B21) 


(B22) 


(B23) 


*( 0 ) = 


(c + v) v - ar 
(c + v) 2 


(B24) 
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appendix c 

REDUCTION OF EQUATION (33) TO A FIRST-ORDER APPROXIMATION 


b 


i 


(c-v)v +ra (c+v)v-ra 

(c-v ) 2 (c+v ) 2 

2 + (c-v)v + ra (c + v) v - ra 
(c - v ) 2 (c + v ) 2 


( 33 ) 


To determine the first-order approximation, the following maximum antici- 
pated values must be assigned to the variables: 

r < 7 x 10 6 meters 
v < 4 x 10 3 meters per second 
a < 5 meters per second squared 
c = 3 x 10 8 meters per second 

When these values are substituted into the denominator, it is obvious that: 


2 >> (c-v)v+ra _ (c-v)v-ra 

(c - v ) 2 (c + v ) 2 


(Cl) 


Therefore, 


b. s -2L 
1 2 


(c - v) v + ra (c+v)v-ra 


(c-v ) 2 


(c +v ) 2 


(C2) 


w f(c-v)v + ra (c+v ) 2 (c+v)v-ra (c-v ) 2 

^ [_ (c-v ) 2 (c+v ) 2 (c+v ) 2 (c-v) 2 _ 


/no\ 
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But 


Therefore, 


where 


and 


Therefore, 


~ w [(cv - v 2 + ra) (c 2 + 2cv + v 2 ) 

2 (c 2 - v 2 ) 2 


- (cv + v 2 - ra) (c 2 - 2cv + v 2 )]. 


(C4) 


c 2 » v 2 . 


(C5) 


_ ~ W [ c 3 v + 2c 2 v 2 + cv 3 - c 2 v 2 — 2cv 3 — V 4 
1 2c 4 


+ rac 2 + 2crav + v 2 r a - c 3 v + 2c 2 v 2 - cv 3 
- c 2 v 2 + 2c v 3 -v 4 + c 2 ra - 2cvra + rav 2 ] , 

= " w _ [2c 2 v 2 + 2c 2 ra - 2v 4 + 2v 2 ra], ( C7 ) 

2 c 4 

= ~ w _ v 2 + ra - — (v 4 - v 2 ra) , 
c 2 c 2 

v 2 + r a » Z— » (C9) 

n 2 


v 2 + ra »-X_ ra- (CIO) 

c 2 


b % (v 2 + ra). (Cll) 

c 2 
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APPENDIX D 

SERVOMECHANISM OF THE TEST SYSTEM 


As discussed in the experimental test system section, it was impossible to 
measure the quality of the Doppler correction scheme in the satellite. The 
signal was therefore looped through the satellite twice, placing a virtual satellite 
(Point A, Figure D-l) on the ground. The only additional circuit functions required 
were two frequency translators to prevent interference of the two transmitted 
signals and of the two received signals, respectively. Appendix A describes the 
operation of the basic servomechanism. 

In the test system, f REF _ is established at the virtual satellite. Figure D-2 is 
a functional block diagram of the error measurement circuitry. Both signals 
were passed through hard-limiters, converting the sinusoidal waveforms to 
pulse waveforms. The positive-going edge of f yg started the time interval 
counter; the positive -going edge of f REF stopped the time interval counter. A 
printout recorded the time interval counts at a rate of six samples per second. 

The initial time interval represents an initial relative phase difference. Any 
subsequent change in interval represents a change in relative phase, or frequency 
difference. 
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Figure D-l— Servomechanism of Test System 
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Figure D-2- Instrumentation of Test System 


TIME 

INTERVAL 

RECORDER 
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APPENDIX E 

REDUCTION OF EQUATION (49) TO A FIRST-ORDER APPROXIMATION 


b 


1 


(c-v)v ira (c + v) v - r a 

(c - v) 2 (c + v) 2 

l + (c-v)v + ra (c + v)v - ra 
(c - v) 2 (c + v) 2 


(49) 


To determine the first-order approximation, the following maximum antici- 
pated values must be assigned to the variables: 

r < 7 x 10 6 meters 
v < 4 x 10 3 meters per second 
a < 5 meters per second squared 
c = 3 x 10 8 meters per second 

When these values are substituted into the denominator, it is again obvious that: 


Therefore, 


1 » ( c ~ v ) v + ra 
(c -v) 2 


(c + v) v - ra 
(c + v) 2 


bj £ -w 


(c - v) v + ra 
_ ( c -v) 2 


(c + v) v - r a 
(c+v ) 2 


(El) 


(E2) 


A comparison of equation (C2) and equation (E2) indicates that b x in the test 
system is exactly twice b t in the operational system, equation (Cll). 

Therefore, in the test system: 


, ^ -2w 2 . 

b £ (v 2 + ra). 

c 2 


(E3) 
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APPENDIX F 


REDUCTION OF A SCATTER DIAGRAM TO A QUADRATIC CURVE 
BY THE CRITERIA OF LEAST-MEAN-SQUARED ERRORS 

Assume that a set of observations, <£,, <£ 2 ,...,<£ , which belong to an un- 
known function, 4> = f (t), at some prescribed points of the independent variable 

t = tj, t 2 t n (FI) 

can be fitted by a polynomial of order n: 

4> = b Q + bj t + . . . + b n t n . (F2) 

The order, n , of the polynomial is chosen a priori; the coefficients, b 0 , b x , b a , 
. . . , b n , are determined by the measurements. 

The unknown coefficients, b.’s, of the polynomial are determined as follows. 
At each point of observation, i , form a "residual" R, so that: 


R. = b. + b, t. + 


+ b n t i n “^i» 


where 4> i is the observed value at t.. 

Take the sum of the square of all residuals: 


(F3) 


Q = 



(F4) 


Q = 7] (b Q + b a t. + . . . + b n t 4>.y, (F5) 

i = 1 
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The quantity ,Q , is by nature positive or in the limit zero. The zero value is 
possible only if all of the residuals vanish independently. This implies that the 
measurements fit the n th order polynomial exactly. This is not the usual case. 
There is, however, a unique value of the bj's for which the sum, Q, becomes a 
minimum. The polynomial associated with these b. 's is the "best fit" of the 
measurements . 

This minimization principle has a unique solution in the form of a linear 
set of equations with a nonvanishing determinant. According to the principles of 
calculus, the condition of minimization requires that the partial derivative 
of Q with respect to any and all b^s vanishes. Repeating equation (F5), 


Q = 2>o+ b > t i+ ... + b n t^ i)2 . (F5) 

i*l 


Taking the derivative with respect to b. yields: 


<b„ + b, t, + ... + b „ (F6) 

j i =1 


Expanding, 


m m m 

<t>i t i i = b o Yl tii + hl Zj tji+1 + ‘ * ’ + b " 

i«l i.l i-l 


where 

j = 0, 1 , 2, . . . , n. 
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In this case n = 2. Equation (F7) can now be expanded as follows for j 
and 2, respectively: 


m m m 

=b 0 m+b 1 2^t.+b 2 ti 2 ' 

i = 1 i = 1 i = 1 


m 


T A*. 




t. 

l 


+ b. 




m m 



Because all measurements of <£. were made at regular intervals of t., 

t. = i At. 

1 

The three preceding equations may now be rewritten as: 


2 > 


= b 0 m + bj A t 


i + b 2 At 2 

i = l 



At 


22 ' i *i =b o at 

iel 


L 


i + b x A t 2 


22 12 + b 2 At3 > 
isl 


m 



i x 1 


m m m m 

A t 2 J i 2 <j> i - b 0 A t 2 2 i 2 + bj A t 3 i 3 + b 2 A t 4 ^ ' i 4 . 

isl i =1 isl 


= 0 , 1 , 


(F8) 


(F9) 


(F10) 


(Fll) 


(F12) 


(F13) 


(F14) 
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These equations, called the "normal equations" of the least-square problem, 
belong to a group of linear equations which have not only a symmetric but also 
a recurrent matrix, characterized by the property 


A jk - A j-i, k+r 


(F15) 


In the preceding set of equations, the b.'s are the unknowns; all other terms 
are constants. The constants involve simple summations which were computed 
on a digital computer. After the constants were determined, a matrix-solving 
routine was used to evaluate the values of b Q , b x , and b 2 . 

Table 1 lists the normalized measured errors, E NM , for the two experimental 
runs. 
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APPENDIX G 

CALCULATION OF ERRORS 


The theoretical errors of the test system, E' NT , are calculated below. The 
errors are determined by using equation (51) and the initial condition listed in 
the experimental section. 


El 


- 2 (v 2 + r a) 


NT 


(51) 


Experiment 1 



9 (4.2) 2 x 10 6 - 6.8 x 10 6 x .25 
3 2 x 10 16 


= -3.54 x IQ” 10 . 


Experiment 2 


E'. 


NT 


9 (2.0) 2 X 10 6 + 3.3 X 10 6 X 4.3 
3 2 x 10 16 


= -7.16 x 10“ 10 . 
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